Abstract. We report the first detection of doubly- 
deuterated methanol (CHD2OH), as well as firm detec- 
tions of the two singly-deuterated isotopomers of methanol 
(CH2DOH and CH3OD), towards the solar-type proto- 
star IRAS16293— 2422. From the present multifrequency 
observations, we derive the following abundance ratios: 
[CHD2 0H]/[CH30H] = 0.2 ± 0.1, [CH2DOH]/[CH30H] 
= 0.9 ± 0.3, [CHgODj/pHaOH] = 0.04 ± 0.02. The to- 
tal abundance of the deuterated forms of methanol is 
greater than that of its normal hydrogcnated counter- 
part in the circumstellar material of IRAS16293— 2422, 
- - - a circumstance not previously encountered. Formalde- 
[ hyde, which is thought to be the chemical precursor of 

• methanol, possesses a much lower fraction of deuterated 
I isotopomers (~ 20%) with respect to the main isotopic 

^ ' form in IRAS16293-2422. The observed fractionation of 
^ \ methanol and formaldehyde provides a severe challenge to 

• both gas-phase and grain-surface models of deuteration. 
Two examples of the latter model are roughly in agree- 
ment with our observations of CHD2OH and CH2DOH if 

I ' the accreting gas has a large (0.2-0.3) atomic D/H ratio. 
^ However, no gas-phase model predicts such a high atomic 
D/H ratio, and hence some key ingredient seems to be 
missmg. 
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1. Introduction 

In the last few years, the study of doubly-deuterated mo- 
lecules in the interstellar medium has gained consider- 
able attention. This field was boosted by the discovery 
of an extremely large amount (D2CO/H2CO ~ 10%) of 
doubly-deuterated formaldehyde in the low mass protostar 



IRAS16293-2422 (hereafter IRAS16293 ; Ceccarelli et al 



1998), a fractionation about 25 times larger than in Orion 
(Turner 1990). Follow-up observational studies of this first 
discovery confirmed this very large degree of deuteration 
in IRAS16293 ( [Loinard etal] 2000), and allo wed a study 
of its spatial distribution (Ceccarelli et al. 2001). Sub- 
sequently, similarly large amounts of doubly-deuterated 
formaldehyde and ammonia have been observed towards 
another very young protostellar core, 16293E, which lies 
i n the same m olecular cloud (L1689N) as IRAS16293 
( Loinard et al, 2001) a nd in the molecular cloud L1689N 
itself ( Ceccarelli et al. 2002). Finally, preliminary results 
of an ongoing project show that it is likely that all low- 
mass protostars present similarly large abundance ratios 
of doubly deuterated formaldehyde wi th respect to H 2CO, 
whereas high-ma ss protostars do not ( [Loinard et al. 2002, 



CeccarelH et al.| 2002). 

All these observations suggest that such a large deuter- 
ation of formaldehyde is produced during the cold and 
dense pre-collapse phase of low-mass protostars. Highly 
deuterated ices are very likely formed via active grain 
chemistry ( Tielens| 1983), stored on the grain mantles, 
and eventually released into the gas during the collapse 
phase, when the heating of newly-formed protostars evap- 
orates the CO-rich ices ( Ceccarelli et al. 2001a). Methanol 
is also commonly believed to be formed on grain surfaces, 
because gas-phase models cannot account for the large de- 



tected abundances of methanol in hot cores ( Menten et al 
1988). If formaldehyde and methanol are produced on 
grain surfaces by simple successive hydrogenations of CO, 
then the reproduction of the abundance ratios between 
deuterated isotopomers and their normal counterparts is 
a crucial test for the grain-surface theory of deuteration 
(e.g. [Charnley et al.| 1997). 



Although the grain picture seems qualitatively con- 
sistent with all the observations so far available towards 
protostars, the nature of the production of deuterated 
molecular species, whether it occurs completely via ac- 
tive grain chemistry (starting from a high D/H atomic 
ratio derived from gas-phase chemistry) or at least par- 
tially via gas-phase formation, is still largely debated (see 
for example Roberts & Millar 2000b). The debate has not 
been settled conclusively because of the relatively small 
body of available observations and the discovery of rel- 
atively large abundances of doubly-deuterated ammonia 
(NHD2/NH3 ~ 0.001; [Roueff et alj 2000) in the molec- 
ular cloud L134N and triply-deuterated ammonia in the 



low- mass protostar NGC1333-IRAS4 (van der Tak et al 



2002) and in the dark cloud Bl (Lis et al. 2002). The 



observed fractionation of ammonia in L134N can be ac- 
counted for by gas-phase models if a high degree of deple- 
t ion of heavy mate rials o nto the grain mantles is assumed 
( [Roberts fc Millarl 2000 b, iRodgers fc Charnley! 2001). The 
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ND3 observations from [Lis et al.[ (2002) can also be ex- 
plained in the framework of gas-phase chemical models 
if the dissociative recombination of partially deuterated 
ions results in a somewhat higher probability for the ejec- 
tion of hydrogen atoms than for deuterium atoms. We 
wish to emphasize that this debate is not merely aca- 
demic, as it involves our understanding of the chemistry 
of the interstellar medium and of ice formation in general 
and deuteration processes in particular. Many observa- 
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tional studies use deuteration processes, which are sup- 
posedly well-understood, to derive key quantities such as 
the deuterium abundance (e.g. in the Galactic Center ; 



Lubowich et al. 2000) or th e degree of ionization (e.g. in 



protostars ; Williams et al.| 1998). The actual state of our 



comprehension of those processes has therefore a large im- 
pact. 

In this Letter, we report the very first detection of 
a doubly-deuterated isotopomer of methanol (CHD2OH), 
with 15 detected lines, towards the low-mass protostar 
IRAS16293-2422. We also report the detection of the 
two singly-deuterated forms (CH2DOH and CH3OD) of 
methanol towards the same object. We compare the de- 
rived fractionation ratios as well as the formaldehyde frac- 
tionation (Loinard et al. 2000) with predictions based on 
active grain chemistry. 

2. Observations and results 

Using the IRAM 30- meter telescope (Pico Veleta, Spain), 
we detected the 23 CH2DOH lines, 6 CH3OD lines, and 
15 CHD2OH lines reported in Table |l|. The telescope was 
pointed at the coordinates a(2000) = 16'^32'"22.6" and 
5(2000) = -24°28'33.0". The observations were performed 
in November 2001 and May 2002. Four receivers were used 
simultaneously at 3, 2, 1.3 and 1.1 mm with typical sys- 
tem temperatures of about 100, 200, 400 and 800 K respec- 
tively. These receivers were connected to an autocorrelator 
divided in up to eight units depending on the setting. The 
telescope beam width varies between 30" at 83 GHz and 
11" at 226 GHz. All observations were performed using the 
wobbler switching mode with an OFF position 4' from the 
center. The pointing accuracy was monitored regularly on 
strong extragalactic continuum sources and found to be 
better than 3". Our spectra were obtained with integra- 
tion times ranging from 70 to 400 min for some CHD2OH 
lines. 

Examples of observed spectra are shown in Fig. |l|. The 
emission was detected at the Vf^sR of the source, namely 
3.9 km s~^. The lines have FWHM of ~ 5 km s~^, very 
similar to the linewidths observed for the HDCO and 
D2CO lines in IRAS16293 ( [Loinard et al| 2000). The mea- 
sured intensities, linewidths and main-beam temperatures 
are reported in Table |l|. 

3. Derivation of the column densities 

We derived the total column density of the three species 
from the rotational diagrams presented in fig. ^j. The ob- 
served scattering is likely due to non-LTE and opacity 
effects. Note that the level column densities are averaged 
on a 10" beam, i.e. the smallest beam of our observations. 



following the suggestion by van Dishoeck et al. (1995) of 
enhanced methanol emission in the central < 10" region 
of IRAS16293. The derived rotational temperatures and 
overall column densities for the three molecules are re- 
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Fig. 1. Examples of CH2DOH, CH3OD and CHD2OH 
lines. The intensities are reported in main-beam bright- 
ness temperature. 
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Fig. 2. Rotational diagrams for the observed transi- 
tions of CH2DOH (diamonds), CH3OD (triangles), and 
CHD2OH (squares). For the sake of clarity the CH3OD 
and CHD2OH points have been shifted by -2 and -4, re- 
spectively. The dashed line is the result of the fit obtained 
by fixing the rotational temperature at 47.5 K. 



ported in Table ||. The CH3OD rotational temperature is 
significatively lower compared with the other two, very 
probably because the detected lines have lower upper 
level energies. We therefore estimated the CH3OD col- 
umn density assuming the same rotational temperature of 
the two other species, namely 47.5 K and find the value 
(2.8 ± 0.3) X 10^"* cm~^. Using the methanol column den- 



sity derived by van Dishoeck et al. (1995) (3.5x10 cm~ 
with a source size of 10"), we finally obtain the following 
fractionation ratios: 

[CH2DOH]/[CH3 0H] = 0.9 ± 0.3 

[CH30D]/[CH30H] = 0.04 ±0.02 

[CHD2 0H]/[CH30H] = 0.2 ±0.1 
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Table 1. Parameters of the observed CH2DOH, CH3OD 
and CHD2OH transitions. The fluxes were derived using 
gaussian fits, and the uncertainty given is ^^stat^ + '^cal^ 
where cr^f^g^^^ is the statistical error and cFq^qI the calibration 
uncertainty (15%). The frequencies for CHD2OH are from 
Mukhopadhyay (in prep.). The line strengths (given in 
Deb^) used to interpret the data are from J. Pearson (pri 



vate communication) for CH2DOH, from Anderson et al 



(1988) for CH3OD, and from |N/lukhopadhyay| (1998) for 
CHD2OH. 



Table 2. 10"-averaged column density N(Molecule) and 
rotational temperature T^ot of CH2DOH, CH3OD and 
CHD2OH respectively. "Assuming Ti.ot=47.5 K (see text). 



Molecule 


N(Molecule) (cm"^) 


T,ot (K) 


CH2DOH 


(3.0 ±0.6) X 10" 


48 ±3 


CH3OD 


(1.5 ±0.7) X lO^"* 


20 ±4 


CHsOD" 


(2.8 ±0.3) X 10" 


47.5 


CHD2OH 


(6.0 ± 2.2) X 10" 


47 ±7 



Frequency 

GHz 



Transition 



Av 



. i4. Discussion and conclusions 



CH2D0H 



89.2512 
89.2754 
89.4079 
110.1054 
133.8473 
133.8729 
133.8818 
133.8929 
133.8974 
133.9302 
134.0655 
134.1124 
134.1854 
207.7808 
223.0713 
223.1073 
223.1283 
223.1311 
223.1311 
223.1537 
223.1537 
223.3154 
223.4223 



CH3OD 



110.1889 
110.2626 
110.4758 
133.9254 
223.3086 
226.5387 



CHD9OH 



83.1292 
83.2895 
83.3036 
166.234 
166.271 
166.297 
166.298 
166.304 
166.327 
166.435 
207.771 
207.827 
207.864 
207.868 
207.869 



2o,2- 
2o,2- 
2o,2- 

9l,8- 

3o,3- 
3o,3- 

32,2- 
32,2- 
32,1- 
32,1- 

3o,3- 

32,2- 
32,1- 
2l,2- 
52,3- 

5o,5- 

52,4- 
54,1- 
54,2- 
53,3- 
53,2- 
52,3- 
52,4- 



10.1 01 
■lo,i 01 
■lo,i Co 
9o,9 01 
■2o,2 ei 

20.2 01 
22,1 01 
■22,1 01 
22,0 Ol 

■22,0 ei 
■2o,2 eo 

■22,1 ©0 

■22,0 eo 

■3o,3 Go 

■4l,4 01 

4o,4 01 

■42,3 01 

44.0 01 

44.1 01 

43.2 01 
43,1 01 

■42,2 ei 

■42,3 Oo 



0.7 
0.7 
0.8 
3.3 
1.1 
1.1 
0.6 
0.6 
0.6 
0.6 
1.2 
0.7 
0.7 
0.3 
0.6 
1.8 
1.4 
0.6 
0.6 
1.0 
1.0 
1.4 
1.8 



17.2 
13.8 

4.5 
83.4 
18.3 
21.7 
33.6 
27.4 
33.6 
27.4 

8.9 
20.2 
20.2 
15.9 
33.6 
35.1 
40.8 
79.4 
79.4 
60.8 
60.8 
40.8 
33.6 



0.19±0.04 
0.16±0.03 
0.31±0.05 
0.42±0.08 
0.60±0.11 
0.60±0.11 
0.34±0.07 
0.67±0.14 
0.30±0.07 
0.37±0.15 
0.67±0.33 
0.43±0.15 
0.26±0.06 
1.01±0.25 
1.13±0.42 
0.51±0.18 
0.54±0.26 
0.41±0.16 
0.41±0.16 
0.58±0.09 
0.58±0.09 
0.38±0.09 
0.55±0.09 



0.04 
0.06 
0.06 
0.05 



0.07 
0.05 
0.08 
0.07 
0.10 
0.06 
0.06 
0.11 
0.17 
0.15 
0.22 
0.08 
0.08 



0.16 
0.10 



li-lo 
2i-2o 
3i-3o 
I1--I0+ 
5n--4i + 



5o 



-f-40-f 



1.6 
2.7 
3.8 
3.2 
3.4 
3.5 



7.8 
10.8 
15.4 

6.0 
26.8 
22.7 



0.10±0.02 
0.17±0.04 
0.30±0.05 
0.55±0.18 
0.09±0.03 
1.32±0.21 



0.04 
0.03 

0.06 
0.06 
0.23 



2o- 
2o- 
2o- 
4o- 

42- - 
43+- 

43- - 

42-I-- 

4o- 
4o- 
5o- 

52- - 
54-- 

53- - 
53-f - 



■lo 01 
lo Co 
Iq 01 

■3o ci 

■32- Ci 

33- ci 
33- ci 

■32+ Ci 

3o Ol 
■3o Co 
■4o ci 

■42- ci 

■44- ci 

■43- ei 

■43-f 01 



1.4 
1.4 
1.4 
2.8 
2.1 
1.2 
1.2 
2.1 
2.8 
2.8 
3.5 
2.9 
1.3 
2.2 
2.2 



16.98 
4.17 
10.33 
26.69 
35.65 
46.54 
46.54 
35.65 
20.04 
13.89 
33.63 
42.59 
68.41 
53.48 
53.48 



0.07±0.02 
0.09±0.02 
0.07±0.01 
0.39±0.08 
0.22±0.07 
0.10±0.08 
0.10±0.08 
0.22±0.06 
0.17±0.10 
0.48±0.13 
0.44±0.11 
0.19±0.11 
0.43±0.13 
0.17±0.10 
0.17±0.10 



0.02 
0.03 
0.03 

0.04 
0.03 
0.03 
0.06 

0.05 
0.09 
0.07 
0.07 
0.05 
0.05 



4.3±0 
2.6±0. 
4.7± 
7.4±0 



4.3±0 
12.2±2 
3.4±0. 
4.8±2. 
6.0±3 
6.4±2. 
4.1±0 
8.4±2 
6.3±2. 
3.1±1 
2.3±0. 
4.8±2. 
4.8±2 



Our most dramatic result is the detection in 1RAS16293 
6Df a form of doubly-deuterated methanol along with 
"Ihe detection of both possible singly-deuterated iso- 
iopomers of this molecule. Up to now, only a ten- 
tative detection of CH3OD has been reported in a 
dow-mass protostellar source ( van Dishoeck et al. 1995). 

ingly-deuterated met hanol has been d efinitely ob- 
served towards Orion ( Mauersberger et al. 1988, where 
ipH30D/CH30H 



1993, where 



2.2±0. 
5.2±0 



2.3±0 
4.5±1 

8.4±2 
1.4±0 
5.5±0 



3.1±0, 
3.2±0, 
2.5±0, 

5.1±1, 
3.2±2, 
3.2±2, 
3.4±0, 

8.7±2, 
4.4±1, 
2.5±1, 
5.6±1, 
3.4±1, 
3.4±1, 



0.01-0.06 and J acq et al 

£;H2D0H/CH30H - 0.04) and SgB2 ( |Gottheb et al 
J979, with CH3OD/CH3OH - 0.01). Equally strikingly, 
2ffe find the deuterated forms of methanol to possess a 
!^otal abundance greater than the main isotopomer in 
i[RAS16293, even without the contribution of the doubly- 
deuterated isotopomer CH2DOD! To date, no other mo- 
decule has been observed to show such extreme deuterium 
'I'ractionation. 

As discussed in the Introduction, the abundances of 
deuterated methanol and deuterated formaldehyde pro- 
Vide a strong test of models involving active grain chem- 
astry. The basic hypothesis behind these models is that 
I'ormaldehyde and methanol form by the hydrogenation of 
-GO accreted onto the grains via reactions with atomic hy- 
rlrogen dTielcns fc Hageij 1982; |Charnley et"ar| l997, here- 



[CH2DOH]/[CH30D] = 20 ± 14 

[CHD2 0H]/[CH2DOH] = 0.2 ±0.1 . 

If the CH3OD column density derived with Trot=47.5 
K is taken, the ratios involving this species change by a 
factor of two, reflecting the increased abundance of this 
molecule. 



|ifter CTR97). The enhanced deuteration is caused by an 
enhanced (atomic) D/H ratio in the gas during the era 
^f mantle formation ( [TielensI 1983). The hydrogenation 
iind deuteration of CO is predicted to form H2CO first 
|C0 ^ HCO H2CO) and subsequently CH3OH (H2CO 
II3CO CII3OH). With some simplifying assump- 
tions, this leads directly to predictions for steady-state 
j-atios of singly- and doubly-deuterated formaldehyde and 
%iethanol to their normal isotopic forms in terms of the 
relative accretion rates of H and D with respect to CO 
"as" free parameters (CTR97). The relative accretion rate 
of H with respect to CO can be derived from the obser- 
vation of the CO/CH3OH and H2CO/CH3OH abundance 
ratios. The predictions for the fractionation ratios then 
depend only on the relative accretion rate of D with re- 
spect to CO, or equivalently on the D/H atomic abun- 
dance ratio in the accreting gas. Figure shows fraction- 
ation ratios predicted by the CTR97 model. In particu- 
lar, the calculated ratios of singly-deuterated isotopomers 
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to normal species are plotted against the analogous ra- 
tios for doubly-deuterated species, both as functions of 
the D/H atomic ratio in the gas. As seen in the upper 
panel, the CH2DOH and CHD2OH observations are com- 
patible with an atomic D/H ratio of 0.2 in the accret- 
ing gas. The CH3OD abundance falls short of the pre- 
dicted value for a D/H ratio of 0.2. However, the gas 
phase abundance of CH3OD in the hot core can be affected 
by gas phase ion- molecule reactions. Specifically, protona- 
tion of methanol by HJ or HsO"*" followed by dissociative 
electron recombination back to methanol will drive the 
CH3OD/CH3OH ratio to the deuterium fractionation of 
molecules in the warm gas, which is very low (CTR97). 
The timescale for this process is some 3x10* yrs which 



2x10", 

We note that this chemical reshuf- 



is comparable to the lifetime of IRAS16293 ( 
CeccarelH eraL| 2000) 



fling of the deuterium will not affect the deuterium frac- 
tionation on the methyl group (eg., CH2DOH, CHD2OH). 
The CTR97 model, however, has some difficulties in ex- 
plaining the observed formaldehyde fractionation ratios 
towards IRAS16293 (Loinard ct al. 2000). In particular, 



the value of atomic D/H compatible with CH2DOH and 
CHD2OH is reasonably compatible with HDCO but re- 
sults in too low a fractionation ratio for D2CO by a factor 
of 5 or so. Possible gas-phase alterations have not been 
considered. The formaldehyde discrepancy suggests that 
at least some of the assumptions in the CTR97 model may 
be too drastic. For example, the CTR97 model contains 
the approximation that only the accreted H, D and CO are 
important in the formation of formaldehyde and methanol 
on the grain surfaces, and that no other reactions compete 
with formaldehyde and methanol formation. 

Very recently, [Caselli et al] (2002) [hereafter CSS02] 
proposed a somewhat more detailed but related model 
for the formation of formaldehyde, methanol, all of their 
deuterated isotopomers, and selected other species on 
grains. In their model, CSS02 consider accretion onto 
grains of H, D, CO and O, followed by a comprehensive 
set of reactions to form H2O, H2, CO2, H2CO, CH3OH 
and all singly- and multiply-deuterated isotopomers of 
these species. Their model differs slightly from that of 
CTR97 in that it specifically includes small differences 
in the barriers to reaction between non-deuterated and 
deuterated species. Yet, in the so-called accretion limit 
and with the same assumption of rapid diffusion rates, the 
CSS02 model should yield approximately the same results 
as the CTR97 model given the same set of chemical re- 
actions and physical conditions. Predictions for methanol 
and formaldehyde fractionation ratios are indeed similar 
to the CTR97 model and substantially the same discrep- 
ancies remain. In particular, for D/H = 0.3, a tempera- 
ture of 10 K, and so-called high-density conditions, CSS02 
agree approximately with our observed fractionation ra- 
tios for CH2DOH and CHD2OH, but produce approxi- 
mately 5 times too much CH3OD, a value similar to that 
of CTR97 with D/H = 0.2 and no subsequent gas-phase 
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Fig. 3. Model predictions (solid lines, adapted from 
CTR97) of abundance ratios between singly-deuterated 
and doubly-deuterated isotopomers with respect to nor- 
mal species are plotted vs one another for methanol 
and formaldehyde. The predictions are obtained as 
functions of the gas -phase atomic D/H ratio. The 
observed ratios of Loinard et al. (2000) are also 



shown. Upper panel: CH2DOH/CH3OH (diamond) and 
CH3OD/CH3OH (triangle) versus CHD2OH/CH3OH. 
Lower panel : HDCO/H2CO versus D2CO/H2CO. 



chemistry. For formaldehyde, CSS02 obtain a fractiona- 
tion ratio for HDCO that is twice the observed value and 
a fractionation ratio for D2CO that is a factor of two below 
the observed value. 

In summary, neither model is in good agreement with 
all of our data. More importantly perhaps, the CTR97 and 
CSS02 models require a D/H atomic ratio in the range 0.2- 
0.3, which is a significantly larger value than can be pro- 
duced by current gas-phase models, even in the presence 
of a large CO depletion (e.g. [Roberts fc Millai] 2000a). Fu- 
ture progress will probably require more detailed chemical 
models in which gas-phase and surface chemistry occur si- 
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multaneously. Although such models arc currently in ex- 
istence, they do not yet contain fractionation processes. 
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